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1 Introduction 

The motion of an electromagnetic wave, through a classically-forbidden region, 
has recently attracted renewed interest because of its implication with regard to 
the theoretical and experimental problems of superluminality. From an exper- 
imental point of view, many papers provide an evidence of superluminality in 
different physical systems]^, |[ |[ ^, ||, |j. Theoretically, the problem of a pas- 
sage through a forbidden gap has been treated by considering plane waves at 
oblique incidence into a plane parallel layer of a medium with a refractive index 
smaller than the index of the surrounding medium, and also confined (Gaus- 
sian) beams, still at oblique incidence]^, ^, ^, |lOj. In the present paper the case 
of a Bessel beam is examined, at normal incidence into the layer (Sees. II and 
III), in the scalar approximation (Sec. IV) and by developing also a vectorial 
treatment (Sec. V). Conclusions are reported in Sic. VI. 

2 The Bessel beam 

An interesting solution of the wave equation is represented by the Bessel beam 
with axial symmetry having, as known the following expression: 

u{p, Tp, z) — AJo{konp sihBq) ex.p{ikonz cos 9q) , (1) 

where A is an amplitude factor, Jq denotes the zero-order Bessel function of first 
kind, p, tp, z are cylindrical coordinates (Fig. la), 6*0 is the parameter of the 
beam (Axicon angle), n is the refractive index of the medium where the beam 
propagates, and /cq is the wavenumber in the vacuum. The beam is independent 
of the angular coordinate ip. The time factor exp(— icjt) is omitted in Eq. ([l]). 

The unusual features of a Bessel beam are that its phase propagates (in the 
z direction) with a velocity v — c/{n cos Oq) larger than c/7i|l^, |l^, and that it 
does not changes its shape during propagation (the amplitude is independent of 
z). The situation is similar to what occurs when only two plane waves interfere, 
the only difference being that the two-wave interference pattern occupies the 
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Figure 1: The coordinate systems 



whole space, while the field is practically limited to a restricted zone of 
space {konpsindo < 2.4, 2.4 being the first zero of Jq). In this connection, it 
is worth noting that the field of Eq. (0) is not properly a beam, since it is not 
limited by a caustic surface, inasmuch as Jq oscillates when its argument tends 
to infinity. 

The meaning of the parameter 6*0 is the following. Let us refer the space to 
a system of Cartesian coordinates x, y, z, (unit vectors i, j, k), such that 

X = p cos ip 
y = p sin ij) 

z = z (2) 

Let us also consider a system of spherical coordinates r, 9, Lp, with origin in the 
origin of the Cartesian coordinates, and the 6 = semiaxis coinciding with the 
positive 2;-axis. 
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Let us now consider a set of plane waves, with the same ampUtude Ad(p, and 
with directions of propagation s = ai+ (3j + 7k making the same angle 6*0 with 
the z-axis. Each wave can be written as 

Aex.p[ikQn{ax + Py + jz)]d(p , (3) 
where the well-known relations hold (see Fig. lb): 

a = sin 9o cos (p , j3 = sin Oq sin tp , 7 = cos 6*0 . (4) 

If we integrate expression (|3|) over df between and 27r, and recall the properties 
of the Bessel function Jq [Q , we arrive at Eq. (|l|) 

l'2-K 

A I exp[ifcon(ax + /3?; + 72)] dip 
Jo 

= A exp{ikonz cos 6q) / exp[iA;o"-p(cos(/?cos + sin(^sin'(/') sin^o]'^^' 
= A 2TTJo{konpsm9o) eKp(ik()nzcos9o) (5) 

If beam (|l|) impinges normally into a plane parallel layer of refractive index n', 
limited by the planes z = and z — d (Fig. 2), the plane waves (^) form an inci- 
dence angle i = 9o which may be larger than the limit angle iq = arcsin(n/n'). In 
this case, the plane waves undergo total reflection, and the layer is a classically- 
forbidden region for the Bessel beam, in spite of the fact that the beam impinges 
normally into the layer. 



3 Total reflection for a Bessel beam 

The preceding expansion of a Bessel beam in plane waves whose directions of 
propagation cover a conical surface of semiaperture 0o immediately indicates 
what happens when the Bessel beam impinges (at z = 0) into a plane surface 
separating two media of different refractive indexes n and n'. Each plane wave 
forms the same incidence angle 9q and, hence, also the same refraction angle 9' , 
satisfying 



n' sin 9' = n sin 9q . (6) 

The refracted waves have the same amplitude, therefore, their superposition u 
is a Bessel beam with parameter 9': 



u{p,ip,z) = A'Jo{kon'psm9')exp{ikon'zcos9') (7) 

Here, it may be interesting to note that, if cos 6*' is real, the phase of the Bessel 
beam for z > propagates in the z-direction, normally therefore to the bound- 
ary. The phase velocity turns out to be c/(n' cos 6*'), i.e. larger than the light 
velocity in the n' medium. If cos 9' is purely imaginary (that is, for larger 
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Figure 2: The layer and the impinging Bessel beam characterised by the Axicon 
angle 9o- 



than iq), the single plane waves composing the incident Bessel beam are in total 
reflection, and therefore give rise to plane refracted waves the phase of which 
propagates parallel to the boundary. However, for the refracted Bessel beam, 
there is no phase propagation for z > 0: there is a sort of stationary field. Here 
as follows, only the latter case will be considered, namely: 

sine' = — sin 6*0 > 1 (8) 
n' 

with the notation 

cos6'' = ir, (r>o). (9) 



4 The tunneling effect in the scalar approxima- 
tion 

If the beam of Eq. (|l|) impinges normally into the layer of Fig. 2, it gives rise: 

i - on the left of the first boundary (z < 0), to a reflected Bessel beam which 
propagates in the negative direction of the z-axis; 

ii - inside the layer (0 < z < d), to two Bessel beams, a "progressive" one and 
a "regressive" one; 
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in - on the right of the second boundary (z > d), to one transmitted Bessel 
beam, which propagates in the positive z-direction. 

The continuity conditions of the total field at the boundaries z = and z — d 
may be easily satisfied with a suitable choice of the complex amplitudes of the 
beams, and of the parameter 9' in the argument of the Bessel beams inside the 
layer. Using standard procedure, let us write the incident field in the form 

u^Pj^Pjz) — A'^2TTJo{konpsm9o)exp{ikonzcos9o) (2 < 0) (10) 
the reflected field in the form 

u'^{p,^p,z) — A"" 2nJo{konpsm9o)exp{—ikonzcos9o) < 0) (11) 

the transmitted field in the form 

u*{p,^, z) — A* 2nJo{konp sin 9o) exp{ikon{z — d) cos 9q) {z > d) (12) 
and, lastly, the progressive and regressive fields, respectively, in the forms 

u'^{p,ip, z) — A^ 2TrJo{kon' psm9') e'xp{—kon' zT) 

u~{p,ip,z) = 27rJo(fco7T.'/9sin6'') exp(fcon'zr) (0 < 2; < d). (13) 

Due to the exponential dependence of u'^ and on z, the progressive and 
regressive beams inside the layer may be denoted as "evanescent" Bessel beams. 

The ratio R = /A^ is the refiection coefficient of the layer for the Bessel 
beam; the ratio T = A*- /A^ is the transmission coefficient. 

By denoting the total field at any point of the space by the boundary 
conditions can be written as 



0- 



'(z = 0+) 



M*°*(Z 

i'°\z = d - Q) = u'°\z = d + Q) 



(14) 



and 



dz 



M*°*(z = 0-) = — u*°*(z = 0+) 



d 



u'°\z = d-0) = -^^ u'°\z = d + 0) 



(15) 



Upon the introduction of Eqs. ([T0|) to (|3|) into Eqs. (|lj) and (|15|), the 
following conditions are found: 
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A' + A'' = A+ + A- 

(A* -yl'^)ncos6'o = in'T{A+-A-) 

A+ei+A^e2 = A* 

in' r{A+ei- A- 62) = A^ncosOo, (16) 

where ei = exp(— fegn'rir), 62 = 1/ei. 

The solution of system (|l^) is easily found and we have[] 



^ = ^(n'r-z.cosW 

4r = ^("'r + zncos0o)r 

^ 4mnT cos6>o 

e2(?T.cos0o + inT)^ — ei(ncos6'o — in'T)'^ 

From the point of view of the field iside to the forbidden region, it should be 
noted that A"*" and A~ are complex quantities that depend on and on the 
geometric characteristics of the system. If we denote the argument of T by 
and if we introduce a quantity $' such that 

^ , n cos 00 
tan$' = -—^ 18 
n i 

we can write 



A+ = |A+|exp[j($T - $')] 
A" |A-|exp[i($T + «>')]- (19) 

Accordingly, the internal total field can be written as 
= u+ + = 27r exp(i$T)>/o(fco"-psin^?o) x 

X [I A+ 1 exp(-fcon'zr) exp(-i<I>') + | A" | exp(fcon'zr) exp(i<I>')] ,(20) 

which shows that, inside the forbidden layer, the total field has a phase <E> = 
$T + with r\[z) such that 



tan,(z) . l^-|exp(fcon-zr)-|A+|exp(-fcon-zr) 
^ |A-|exp(fcon'^r) + |A+|cxp(-fcon'zr) 

- tanh[fcofi'r(d - z)\ tan , (21) 



^The evaluation of the reflection coefBcient, which maybe derived from Eq. (p^, is of no 
interest in the present paper. 
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which propagates in the direction of the positive z-axis. 

Equation may be used to evaluate Igradry] — d-q/dz, the wavelength 
A' = A'(z) and the phase velocity v' = v'{z): 



Igradryl 



kon cos 9o cos^[r7(z)] 
cosh^[konT{d- z)] 

2tt _ Ao cosh^[A:on'r(d- z)] 
Igrad?]! n cos 00 cos^[r/(z)] 

Lu c cosh^[A:on'r(d — z)] 



|grad?7| rt cos 6*0 cos2[77(z)] 



(22) 



where Aq = 27r/fco denotes the free space plane- wave wavelength. It is interesting 
to note that for z ~* d, v' c/{n coado), which is equal to the phase velocity 
of the trasmitted (z > d) and incident (z < 0) fields. This is a phenomenon 
similar to the one reported in Ref . fiot] . 

The phase difference A$ of the total internal field at z = from z = is 
given by 

A$ = arctan[tanh(A;on'rrf) tan $'] , (23) 

while, as to the transmitted field, it can be noted (see Eq. ([T^)) that its phase 
at z = d is 



— arctan 



cos^ 00 - , 

tanh(fcor7, i d) 



2nn'T cos 9i 







(24) 



Thus, $T is equal to the phase shift of the transmitted beam at z = d with 
respect to that of the incident field at z = 0. By comparing Eqs. ( ^ ) and (p^, 
it appears that <I>t does not coincide with the phase gained by the internal field 
in passing from z = to z = d. 



5 The vectorial treatment 

In this Section we use the vectorial algorithm for analysing the propagation of 
Bessel beams. 

5.1 Quasi-TE and quasi-TM beams 

For a vectorial treatment of the propagation of a Bessel beam it is sufficient 
to consider, for example, the function of Eq. (^ as the tangential component 
of the electric field E. Then, the Maxwell equations allow us to determine the 
longitudinal component of E, and the magnetic field H as well. The longitudinal 
component of E turns out to vanish on the z-axis, for p — 0, where the tangential 
component has its maximum. Thus, such field (E, H) may be named quasi-TFi. 
If the field (E, H) impinges normally on the layer of Fig. 2, it gives rise, as in 
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the scalar approximation, to a reflected Bessel beam for z < 0, to a transmitted 
beam for z > d, and to a progressive and a regressive Bessel beams for < z < d. 
The progressive and regressive beams are evanescent if Eq. (||) holds. 

Analogously, we can consider the function u of Eq. as the tangential 
component of a magnetic field H', then from the Maxwell equations we can 
derive the longitudinal component of H' and the associated electric field E'. 
The field (E', H') may be named quasi-TM. 

The complex amplitudes of all the above beams may be determined by im- 
posing, at the two boundaries, the continuity conditions for the tangential com- 
ponent of both the total electric and magnetic fields. This allows to determine 
in both cases the total field inside the forbidden region, its wavelength and its 
phase velocity. 

Since the treatment is a little cumbersome, we report that different wave- 
length and different phase velocity are found for the quasi- TE and quasi-TM 
cases; here we limit ourselves to develop the simpler analysis in the TE and 
TM cases. To this end, let us consider a plane TE and a plane TM waves with 
a direction of propagation s* = ai -|- /3j -I- 7k = sin 9q cos ip'i + sin 9o sin -|- 
cos6'ok. Let them reflect and refract through the first boundary, then reflect 
and refract at the second boundary. Lastly, we integrate the three Cartesian 
components with respect to between and 27r. 

5.2 The TE case 

This case has already been treated in Ref.|^, but in the particular case of 
if — Q. Here, we have to generalise the results obtained therein. 
For the incident field E* = E^'i +Ey\ (E^ = 0), let us put 

Ex = exp[ifcon(a;sin6'o cosiy9 -I- j/sin6'o sin(^ -I- zcos^o) ]'^'^ 

Ey — Ay eyip[ikon{x sin 9q cos if + y sin 9o sin if + z cos 9o)]dip . (25) 

Since E* is normal to s% the following relation must hold: 

A^a + Ayf3 = . 

Hence, 

Ax cos ip + AySinip = 0. (26) 

Equation ( |26| ) shows that Ax and Ay depend on ip. The solution of Eq. (p6|), 
which remains finite for any value of if, is 

Ax = A sin if 

Ay = -Acosip, (27) 

where A is a constant. 
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Figure 3: Amplitude (a) and phase (b) of the transmission coefficients T in 
the TE case (solid lines) and in the TM case (dotted lines) plotted vs d for two 
values of the parameter 6*0- Other parameter values are: u = 60, n = 1.5, n' — 1. 

At this point, insertion of Eqs. ( p7| ) into Eqs. ( p5| ) and integration with 
respect to ip yields (by putting a = 271 A) 

Ex,beam = i tt siiitp Ji{konp sihOq) e:>q) {ikouz COS 9o) 

Eyfieam ~ ~i a cosiljJi{konp sui 0q) e'xp {ikouz COS 9o) , (28) 

where Ji is the Bessel function of the ffi'st order. 

The magnetic field H' — Hx'i + Hyj + i/^k related to the E' field of the 
plane wave is found to be 

H,j. = — cos 9q cos (p exp[ikQn{x sin 9q cos f + y sin 9q sin (p + z cos 9q) ]dip 
Z 

Hy — — COS 9 Q sin if eyi\)[ikQn{x sin 9 Q COS (p + y s\n 9 Q sin Lp + z COS Of)) ]dLp 
Z 

Hz = — — sin exp[ ifcon(a; sin^o cos + y sin ^0 sin + z cos ^o) , (29) 
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where Z — Zq/u^ and Z^ is the free space impedance. Hence, upon integration 
over we obtain 



% Qj 

Hx,beam — COS 6q COS ipJi [k^np sin ^o) exp {ikonz cos ^o) 
1 d 

Hy^heam = — COS 9q sin ipJi {koTip sin ^o) cxp {ikonz cos ^o) 

HzMam = - ^sinOo Jo{kQnpsmdo)exp {ikonz cos9o) . (30) 

A field described by Eqs. (|2^) and ( pO| ) can be denoted as a Bessel beam (of 
tlie first order), and lias properties very similar to those of the Bessel beam of 
Eq. (^). However, only the longitudinal component of the magnetic field is of 
the type of Eq. (|^), that is, picked on the axis at p = 0; all other components 
for p = vanish. A TE plane wave like the one described above gives rise, 
at the incidence on the layer of Fig. 2, to a reflected first-order Bessel beam 
and to two transmitted "evanescent" Bessel beams. One of these is progressive 
and the other regressive, with coefficients A'^ and A~ given by the first two 
Eqs. (0). At the second boundary of the layer, a transmitted first-order Bessel 
beam originates, with amplitude A* given by the third Eq. (pT|). Accordingly, 
the wavelength Xte and the phase velocity vte inside the layer are again given 
by Eqs. (H). 



5.3 The TM case 

The TM case may be treated in a way analogous to that of the TE case, arriving 
at similar results. The only difference lies in the fact that n must be replaced 
by n', and vice versa, in the coefficients of system (^) (not in the propagation 
factors ei and 62), analogously to what happens in the Fresnel formulas for the 
reflection and transmission coefficients of a real plane wave at a plane interface. 
Accordingly, the expression of , A^ , A* and must be replaced by 



A+' 




_ ^2 


~A^, 


^ TM 


2nT 


A-' 




_ ei 




^ TM 


2nT 










^ TM 


— Ttm 



tan 



(nF — in' cos 6*0) 



Ainn'T cos 9n 



62 (ri' COS0Q + mF)2 — ei(n' cos 0o — inV)^ 



n cos 

7lF 



(31) 



6 Conclusions 

On the basis of the vectorial treatment we are now able to compare the TE and 
TM cases. In particular we can conclude that: 
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Figure 4: Normalised velocities vte/c (solid line) and vtm/c (dotted line), as 
a function of d, for the same parameters values as in Fig. 3 and for 6*0 = 45°. 

1 - The phase shift A<1) of the transmitted TM beam (Eq. (|3|) with $" instead 

of $') at z = d with respect to the incident beam at z = is different 
from that of the TE case. 

2 - The transmission coefficient T for the TM case is different in amplitude 

(Fig. 3a) and phase (Fig. 3b) with respect to that of the TE case. Con- 
sequently, an incident Bessel beam formed by a TE component and by a 
TM component gives rise to a transmitted Bessel beam with a different 
polarisation. 

We note that the amplitudes are slowly varying functions of 6*0, for both 
TE and TM cases. On the contrary, the phase in the case of TE differs 
greatly with respect to the TM case, and both of them greatly vary with 

3 - The phase velocity vtm and the wavelength \tm inside the layer in the 

TM case are different from those of the TE case. 



nAo cosh [kQn'T{d 



n'^ cos 6*0 cos2[77tm(2:)] 

nc cosYv' [kQin' T {d — z)] 
n'^ cos 6*0 cos2[?7TAf(^)] 



Vtm = 2 , , , /I „„„2r„ 731 



where rjTM{z) is such that tan[77TA./ (z)] = — tanh[A;on'r(rf — z)]tan<i>". 
Figure 4 shows the normalised phase velocities vtm/c (which is equal to 
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Xtm/^o) a-iid vte/c (equal to Ate/Aq). Wc note that, immediately after 
the first boundary the motion is extremely fast since the effect of the 
second boundary, which originates the regressive (or "anti-evanescent") 
beam, is negligible. Wc wish to recall that, in the absence of the second 
boundary, the phase velocity is infinite. 
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